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Abstract. We consider a spatially homogeneous and isotropic cosmological model 
where Dirac spinors are coupled to classical gravity. For the Dirac spinors we choose 
a Hartree-Fock ansatz where all one-particle wave functions are coherent and have 
the same momentum. If the scale function is large, the universe behaves like the 
classical Friedmann dust solution. If however the scale function is small, quantum 
effects lead to oscillations of the energy-momentum tensor. It is shown numerically 
and proven analytically that these quantum oscillations can prevent the formation 
of a big bang or big crunch singularity. The energy conditions are analyzed. We 
prove the existence of time-periodic solutions which go through an infinite number 
of expansion and contraction cycles. 
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1. Introduction 



The nonlinear coupling of gravity as described by Einstein's equations to matter 
formed of quantum mechanical waves has led to surprising physical effects and has in- 
cited the development of interesting mathematical methods needed for their analysis. 
T.D. Lee et al [12] studied a complex massive scalar field and found spherically sym- 
metric, static soliton solutions, interpreted physically as "mini-soliton stars." In the 
seminal papers [H0 E], D. Christodoulou analyzed the spherically symmetric gravita- 
tional collapse of a real massless scalar field. The study of more fundamental quantum 
mechanical matter formed of Dirac wave functions was initiated in [10] , where spheri- 
cally symmetric, static soliton solutions of the coupled Einstein-Dirac (ED) equations 
were constructed numerically. Moreover, it was shown that the Einstein-Dirac-Maxwell 
and Einstein-Dirac- Yang/Mills equations do not admit spherically symmetric, static 
black- hole solutions (see [11[ Part II] and the references therein as well as [H [2]). 
The resulting physical picture for the gravitational collapse is that after a horizon has 
formed, the Dirac particles necessarily either disappear in the black hole or escape 
to infinity, but they cannot remain in a bounded region outside the horizon. More 
recently, in the independent works (3j [TBI [15] existence results were proved for the 
spherically symmetric, static ED solitons, using quite different mathematical meth- 
ods. 

For time-dependent, spherically symmetric ED systems, the numerical results of [T7] 
as well as the improved and more detailed numerical analysis in [2Q|. [19] show that 
a "cloud" of Dirac particles either dissipates to infinity or else a black hole forms, 
in which case the Dirac wave partially disappears inside the horizon and escapes to 
infinity. Unfortunately, apart from the dynamical stability analysis in [10] . there are 
no rigorous results on time-dependent, spherically symmetric ED systems. The main 
problem is that, compared to the massless scalar field in [4], the rest mass and the fact 
that the Dirac spinors have several components cause additional difficulties. This was 
our main motivation for searching for a simpler time-dependent ED system which is 
easier to analyze but nevertheless reveals new effects on the nonlinear dynamics. This 
led us to considering spatially homogeneous and isotropic systems. By separating 
the spatial dependence, one can reduce the system to nonlinear ordinary differential 
equations in time. In the recent paper [8], it was shown numerically that in such 
systems, quantum oscillations of the Dirac wave functions can prevent the formation of 
a big bang or big crunch singularity. In the present work, we give a detailed derivation 
of the underlying equations and present a rigorous analysis of all quantum effects. 

More specifically, we consider the simplest possible cosmological model where the 
metric is given by the homogeneous and isotropic line element 



where dQ, 2 is the line element on S 2 , and R(t) is the so-called scale function. The 
parameter k can take the values —1, or 1, corresponding to an open, a flat and a 
closed universe, respectively. We will mainly focus on the closed case, but we will give 
a brief outlook on the other cases as well. As matter we consider Dirac particles of 
mass m, as described by solutions of the Dirac equation (D — m)\I/ = 0. Here the 
Dirac operator clearly depends on the metric. On the other hand, the spinors enter 
the Einstein equations via the energy momentum tensor of the wave functions. This 
mutual coupling gives rise to a system of nonlinear partial differential equations. 
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In order to keep the equations as simple as possible, we want the whole system to 
be again homogeneous and isotropic. This forces us to consider several Dirac particles, 
as we now explain in the closed case (the argument in the other cases is similar). Then 
the hypersurfaces t = const in (jl.ip are isometric to three-dimensional spheres. This 
makes it possible to separate the spatial dependence of \f by the ansatz 

with complex functions a and /3, where ip\ is an eigenspinor of the Dirac operator T> s s 
on S 3 corresponding to the eigenvalue A. The Dirac operator V S 3 has the eigenval- 
ues A = ±|,±|, . . . with multiplicities N = A 2 — 1/4. As the symmetry group SO(4) 
of S 3 acts transitively on the A-eigenspace, we can retain the SO (4)-symmetry only 
by taking an anti-symmetric product of all the wave functions which form a basis of 
this eigenspace. Thus, similar to the fermion configuration in an inert gas, we fill the 
whole "shell" of states of momentum A. The number of particles of our universe equals 
the dimension N = A 2 — 1/4 of the eigenspace. 

We point out that with the ansatz (|1.2p . all one-particle wave functions have the 
same time dependence, thus forming a coherent macroscopic quantum state. Due to 
the analogy to a BCS state in superconductivity or to the ground state of a Bose- 
Einstein condensate, we can thus imagine the fermionic many-particle state as a spin 
condensate. The coherence of all wave functions seems essential for the prevention of 
space-time singularities. For example, one could consider an alternative ansatz where 
one fills several shells corresponding to different momenta Ai, . . . , Xl- Apparently, this 
ansatz would make the prevention of a singularity less likely. Thus the specific features 
of our model should be attributed to a spin condensation effect, leading to a coherent 
macroscopic fermionic state. 

With the above ansatz, the Einstein-Dirac equations reduce to ODEs involving the 
scale function R(t) and the above complex functions a(t) and /3(t). Combining a 
and /3 to a two-spinor and taking the expectation values of the Pauli matrices, we can 
describe the dynamics of the Dirac wave functions by a vector i;£K 3 . After a suitable 
orthogonal transformation w = Uv, the Einstein-Dirac equations take the simple form 



w = d A w , R 2 + k = 7T a/A 2 + m 2 R 2 w\ 

R z 



(1.3) 



where the vector d is given by 



The ^-equation in (|1.3p describes a rotation of the Bloch vector w around the axis d, 
which is itself moving in time. In particular, one sees that the length of the Bloch 
vector is constant. We choose the normalization convention 

\w\ = A 2 - - = N. (1.5) 
1 1 4 

In the limiting cases A = and m = 0, our equations reduce to the the well-known 
Friedmann equations for the dust universe and the radiation universe, respectively. 
Moreover, for large R the universe behaves classically as in the dust case. However, 
near the big crunch or big bang singularities, quantum effects change the situation 
dramatically. More precisely, the Dirac equation on the left of (|1.3p can cause the 
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vector w to "tilt", with the effect that \wi\ in the right equation of (jl.3p can become 
small. If this happens, R can become zero even for small values of R. As we shall see, 
in this case R changes sign, thus preventing the formation of a big bang or big crunch 
singularity. We refer to this effect as a bouncing of the scale function. Our main result 
is to give a rigorous proof of bouncing and to prove the existence of time-periodic 
solutions: 

Theorem 1.1 (Existence of bouncing solutions). Given A G {±|,±|,...} and 
5 > as well as any radius R > and time T >; there is a continuous three-parameter 
family of solutions (R(t),w(t)) of the system (|1.3p defined on a time interval [0, T] 
with T > T > having the following properties: 

(a) At t = and t = T, the scale function has a local maximum larger than R >; 

R(t) > R> , R(t) = , R(t) < . 

(b) There is a time tb ounce £ (0, T) such that R is strictly monotone on the inter- 
vals [0,tb ounce ] and [t bounce ; T]. Moreover, the scale function becomes small in 
the sense that 

R{tbounce) < 5i?> . 

By increasing the parameters T> and i?>, we can arrange that the bouncing solution 
exists for arbitrarily long times and that the scale function becomes arbitrarily large. 
With the parameter 5, on the other hand, we can make the scale function at the 
bouncing point as small as we like. Thus the solutions behave qualitatively as shown 
in Figure [2] on page [TBI 

Theorem 1.2 (Existence of time-periodic solutions). Given A 6 {±|,±|,...} 
and 5 > as well as any radius R > and time T>, there is a one-parameter family 
of solutions (R(t),w(t)) of the system (jl.3p defined for all t € M with the following 
properties: 

(A) The solution is periodic, i.e. 

R(t + T)= R(t) , w(t + T) = w(t) for alltGR, (1.6) 
and every T > with this property is larger than T> . 

(B) inf R{t) < 5 i?> , i?> < sup R(t) . 

R R 

At first sight, this theorem might seem to be in conflict with the Hawking-Penrose 
singularity theorems, which state that a contracting universe must form a space-time 
singularity. However, there is no contradiction because our time-periodic solutions 
violate the strong energy condition at the minimum of the scale function (for details 
see Section l4~3l) . 

The paper is organized as follows. After introducing our model (Section [2]), in 
Section[3]we give a qualitative analysis including scalings, a numerical example, explicit 
approximate solutions and a discussion of the probability for bouncing to occur. We 
proceed with the existence proofs for bouncing solutions (Section |4]) and for time- 
periodic solutions (Section [5]). In Section [6j we briefly discuss the cases of a flat 
and open universe. The appendices provide preliminary and more technical material 
related to Dirac spinors in curved space-time and the derivation of the homogeneous 
and isotropic Einstein-Dirac equations. 
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2. The Model 
The Einstein-Dirac (ED) equations read 

1 
2 



R) - ^R 5) = 8ttk Tj , {V - m)# = , (2.1) 



where T l - is the energy- momentum tensor of the Dirac particles, n is the gravitational 
constant, T> denotes the Dirac operator, and ^ is the Dirac wave function (we always 
work in natural units h = c = 1). For the metric we take the ansatz of the spatially 
homogeneous and isotropic line element 

ds 2 = dt 2 - R 2 (t) ( + r 2 d$ 2 + r 2 sin 2 dip 2 

V/ ( r ) 

where t 6 R is the time of an observer at infinity, (p) G (0, ir) x [0, 2ir) are the 
angular coordinates, and r is a radial coordinate. The function R(t) is referred to as 
the scale function. In the cases of a closed, open and flat universe, the function / and 
the range of the radial variable r are given by 



closed universe: /(r) = yl — r 2 , r G (0, 1) 

open universe: f(r) = yl + r 2 , r > 

flat universe: /(r) = 1 , r > . 

We combine these formulas by writing 



f[r) = \/l- kr 2 with fee {1,-1,0}. 

The Dirac equation in a homogeneous and isotropic space-time is derived in Appen- 
dix O By separating the spatial dependence in the form (jl.2p . we obtain a coupled 
system of ODEs for the complex- valued functions a and /3, 



% — 



dt \P 



X/R -m U ■ (Z - Z) 



Here the quantum number A describes the momentum of the Dirac particle, whereas ip\ 
is an eigenspinor of the spatial Dirac operator. In the closed universe, A can take 
the quantized values ±|,±|, ... (see [HI Appendix A]), whereas in the flat and open 
universes, A can be any real number. We always normalize the spinors according to 

\a\ 2 + \P\ 2 = \ 2 -\. (2.3) 

Note that $ is not homogeneous and isotropic, because its momentum and spin 
orientation distinguish a spatial direction, and because its probability density |^| 2 is in 
general not constant on the hypersurfaces t = const. In order to obtain a homogeneous 
and isotropic system, we consider a family of wave functions (all of the form (|1.2|) with 
the same A and the same (a(t), /3(t))) and build up a complete fermion shell (for details 
see Appendix [B]) . A direct calculation yields that the non- vanishing components of 
the energy- momentum tensor are (see Appendix IB"!) 



r ° = r- 3 



2A — 
m(\a\ 2 - \p\ 2 ) - — Re(a^) 



o \ _ 

X = Tl = T£ = R- 3 — Re(aP). 



(2.4) 
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Similarly, a short calculation yields for the Einstein tensor Gr{ (see [14} Box 14.5, 
eqns (5)]) 

z^o o ^ nr — nti — — o i ^ 

^0 - d R2 ' W - W - - ^ ^ H ' 

and all other components vanish. Thus the Einstein equations Gri = 8ttkTI reduce to 
the two ordinary differential equations 



R 2 + k 
3 — —5 — = 8ttk R 
R l 



2A — 

m(H 2 - |/3| 2 ) -^Re( a /3) 



(2.5) 



6^ + 6 = 8™ R- 3 m (\a\ 2 - |/3| 2 ) . (2.6) 

R R z ' 

Differentiating the first Einstein equation (|2.5p and using the Dirac equation (|2.2p . 
we find that the second Einstein equation (12. 6p is automatically satisfied. Thus this 
equation can be left out. Furthermore, the scaling 

R—j-AR, t -)■ At , m—tA~ l m, A -> A , (a, /3) ->■ (a, /3) , k -> A 2 k 

allows us to arbitrarily change the gravitational constant k, (while preserving the nor- 
malization (I2.3p ). Thus for convenience we may choose 

(2.7, 



Then the Einstein-Dirac equations become 

^(H 2 -|/3| 2 )-| 



R 2 + k = ^ (|a| 2 - |/3| 2 ) - ^ Re(a/3) . (2.8) 



For the analysis of the system of ODEs (|2,2p and (j2.8j) . it is convenient to regard 
the spinor (a, f3) as a two-level quantum state, and to represent it by a Bloch vector v. 
More precisely, introducing the 3-vectors 

/O'K/OL and ^=^ ei - 2me3 m 

(where a are the Pauli matrices, and ei,e2,e% are the standard basis vectors in M 3 ), 
the ED equations become 

■n 1 

v = bAv, R z + k = -—b-v 

2R 

(where 'A' and '•' denote the cross product and the scalar product in Euclidean M 3 , 
respectively). To further simplify the equations, we introduce a rotation U around the 
e2-axis, such that b becomes parallel to ei, 

2 



Ub = — \A 2 + m 2 R 2 e x . (2.10) 
R 

Then the vector w := Uv satisfies the equations (|1.3p and (II. 41) . We refer to the two 
equations in (11. 3p as the Dirac and Einstein equations in the Bloch representation, 
respectively. Our normalization convention (|2.3p implies that the length of the Bloch 
vector satisfies (11.51). 
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3. Qualitative and Numerical Analysis 

3.1. Limiting Cases and Scalings. Let us discuss the ED equations in the Bloch 
representation (jl.3p and (jl .4|) . The Dirac equation describes a rotation of the Bloch 
vector w around the axis d, which is itself moving. The direction of the rotation axis 
is described by the quotient 

|tfe| Ami? \R\ 

N ~ 2(A 2 + m 2 i? 2 )l ' 

If this quotient is close to zero, the rotation axis is almost parallel to the ei-axis, 
so that w\ is a constant. Conversely, if the quotient is large, then the rotation axis 
is parallel to the e2-axis. This has the effect that w\ has an oscillatory behavior. 
This leads to "quantum oscillations" of the energy-momentum tensor in the Einstein 
equations in (|1 .3j) . In physical terms, these oscillations can be understood as the 
Zitterbewegung of Dirac particles. In the intermediate regime when the quotient (13. II) 
is of the order one, the motion of the Bloch vector resembles the precession of a 
gyroscope. We denote the corresponding length scale by i?tilt 5 



= 1 . (3.2) 



M2J 
\di\ 

In the regime R 3> R qu when w\ is constant, the Einstein equation in (|1.3p decouples 
from the Dirac equation. The right side of the Einstein equation in (jl.3p is a monotone 
decreasing function of R. In the case of a closed universe, this implies that R becomes 
zero at a radius R max determined by 

^ iax - wi =: w limax . (3.3) 



VA 2 + m 2 R 2 ma 

This corresponds to the scale when the expansion area of the universe ends, and a 
contraction begins. In the cases of an expanding open or flat universe, it is still useful 
to introduce -Rmax by (13. 3p . It is then the length scale where matter ceases to be the 
driving force for the dynamics of the universe. 

The term A 2 + m 2 i? 2 gives rise to yet another length scale, which we denote by R qu , 

R q u = ~. (3.4) 
m 

If R ^> R qu and R ^> Rait, the Einstein equation is well approximated by the Fried- 
mann equation for dust, 

R(t) 2 + k = — ^ with c = —mwi . (3-5) 
R(t) 

If on the other hand R <C R qu and it! S> -Rtiit , we obtain the Friedmann equation for 
radiation, 



R(tf + k = — %r with c = -A 2 ^i . (3.6) 

R{ty 



The transition from dust to radiation is described by wave mechanics. This is why R qu 
can be regarded as the scale where quantum effects come into play. 

Recall that our equations involve the two parameters A and m. Moreover, a solution 
is characterized by the initial orientation of the vector w in R 3 , giving rise to two 
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additional parameters. Thus our systems admit a four-parameter family of solutions. 
Here it is most convenient to parametrize the solution space by the four quantities 



i? max , A , m and 6 max = arctan — 



(3.7) 

R — ^max 



w 2 

The parameters R qu and wi, m ax are then determined by (|3.4p and (|3.3p . We are 
interested in the situation when a classical universe gets into the quantum regime (or 
similarly if a quantum universe expands into a classical universe). Therefore, we only 
consider the parameter regime where 

Rma,x ~> Rqu, -Rtilt • (3.8) 

Then (|3.3() simplifies to 

Wl,max ~ • (3.9) 

m 

Let us now follow the solution starting at R = R max through an era of contraction up 
to the radius Run- We can then approximate w\ by the constant Wi >max . Moreover, 
for the computation of Rtm we may neglect the summand +k in (jl.3p . Thus we may 
approximate (|3.ip by 

„ 5 

X 2 m 2 \w hmax \ _ mR mSiX ( ^ + I&\ ~2 ^ ^ 



2(A 2 + m 2 i2 2 )§ 2 A 3 V ^qu 



(h 

where in the last step we used (|3.4p . This expression becomes larger if R decreases, 
meaning that quantum oscillations become stronger. Since we want R tnt to exist, we 
are led to demanding that 

^F>1. (3.11) 
Comparing with (13.10P and f)3.2[) , we conclude that 



R q » < Rm * • (3-12) 

To summarize, we always consider the quantities in (|3.7p as our free parameters, which 
we always choose in agreement with (|3.8p and ()3.1ip . Then the scales R ma ,x, Rtnt 
and i? qu automatically comply with (|3.12p . 

3.2. Construction of Approximate Solutions. The above considerations yield a 
method for constructing explicit approximate solutions. We now describe this method. 
In Section 13.31 we shall illustrate it by a numerical example, while in Section 13.41 we 
will use it to get information on the probability of bouncing. For simplicity, we restrict 
attention to the closed case k = 1, although the methods apply similarly in the open 
and flat cases. In short, we begin at t = with R{0) = R m a.x, adjusting parameters 
such that R(0) = and R(0) < 0. We then approximate the Dirac equation by 
keeping the rotation axis w fixed in the direction of e\ (Era I). Then w\ is constant, 
so that R(t) goes over to a solution of the Friedmann equation for dust (|3.5p . We 
use this approximation up to a time t = ttilt when the scaling function attains the 
value .R(itiit) = -Rtiit> where quantum oscillations become relevant. At the time t = t t iit, 
we instantaneously tilt w into the direction e 2 by omitting the term proportional to e\ 
in (jl.4p . From then on, we again keep the rotation axis w in the fixed direction e 2 
(Era II). This approximation of instantaneous tilt is exact in the limiting cases R 3> 
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R t [\ t and R <C Rtnt- For the validity in the intermediate region R « i? ti i t see the last 
paragraph in Section 13.31 

We again take the quantities in (|3.7p as our free parameters, whereas R qu and li^max 
are determined by (|3.4p and (|3.3p . In Era I, we begin at the maximum point R = 
Rma,x 3> Rtnt, Rqu- Thus in the Einstein equation in (jl.3p we replace wi by the 
constant 

wi ]t nax' Using the approximation (j 3 . 9 j) . we obtain the initial value problem 

R + 1 = -77 wi,max , w = 2meiAw (3.13) 

R 

R(0) = i? max , w(Q) = (wi,max,/ocos(^ max ), / osin(0 max )) , (3.14) 

where according to our normalization convention (|1.5|) we chose 



p = y/N*-v% max . (3.15) 

The -R-equation in (|3.13p is the standard Friedmann equation for dust. It can be solved 
by separation of variables and integration (see for example [13} page 138]). The vector 
w(t) is given explicitly by 

w(t) = (w ljmax ,pcos(2mt + max ),psin(2mt + </> max )) • 

We follow this solution up to the time t = i t nt where i2(ttut) = Rtnt- Note that, 
according to (|3.12p . R always stays larger than R qn , making it unnecessary to consider 
the equations (13. 6p for the radiation dominated universe. 
In the following Era II, we approximate the system by 

1 \ttiR 
R 2 + ± = -jpVX 2 + m 2 R 2 w u ^=- x2 + m 2 R 2 e ^A^ (3.16) 

-R(ttiit) = -Rtilt , w(*tilt) = (^l,max,PCOs(0tilt),PSin(0 ti it)), (3.17) 

where we set 

0tilt = 2m ttilt + 0max • (3.18) 

As the u;-equation in (13. 16|) describes a rotation around the e2-axis, the component W2 
is constant. Therefore, it suffices to consider the (ei, e3)-plane. We describe the 
rotation of the vector (1^1,103) in this plane by an angle 9, i.e. 

'wi(t)\ _ ( cos8 sm6\ /wi(t t iit) N 
y w 3(t)J \-sm9 cos6y \w 3 (t ti i t ) / 

Then the u;-equation can be rewritten as 

" R 



(3.19) 



= arctan 

dt 

having the explicit solution 

6(R) = arctan 



Rqu 



Rtnt 


— arctan 


" R ' 




_ Rqu _ 




_ Rqu _ 



9(t m ) = 0, (3.20) 



(3.21) 



This can now be used in (|3. 16[) to obtain the differential equation 



R 2 + l = -- — 2 (w hmax cos 6 (t) + p sm(^ tilt ) sin 9 (t)j. (3.22) 

Clearly, this nonlinear equation has a local solution with initial condition -R(ttiit) = 
i?tiit, which can be maximally extended until 1/R(t) diverges at the moment of a big 
crunch. 
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Figure 1 . The exact solution and the approximation of immediate tilt 
in the case m = 21.5, A = 3/2, R max = 10, R ti i t 0.160, R qu ps 0.0697, 

Wl,vaax ~ -0.465. 



3.3. Numerical Results. We now consider the approximation of instantaneous tilt 
in a typical example. On the left of Figure [H a numerical solution of the ED equa- 
tions (jl.3p and (|l,4p is plotted, where the universe begins in an era of classical con- 
traction. At t ~ 15, quantum oscillations come into play, which trigger a bouncing 
of the scale function leading to a new era of classical expansion. The function w\ is 
almost constant in the classical eras, but is oscillatory in the quantum regime. In the 
approximation of instantaneous tilt, w\ is exactly constant up to the time itnt an d 
oscillates afterwards. As one sees in Figure [lj this is a quite good approximation (in 
the upper left plot for R(t), the approximate solution looks just like the exact solution 
until after the bouncing). 

Clearly, the limitation of the approximation of instantaneous tilt is that it does 
not take into account the "precession" of the Bloch vector around the moving axis of 
rotation w in the region R ~ i?tiit- Generally speaking, the approximation is better 
if the expression on the left side of (13. lip is larger (in the example of Figure HJ this 
expression has the value ~ 63.7). A precise error estimate goes beyond the scope of 
this paper, because it sees difficult to describe the above precession quantitatively. 

3.4. The Probability of Bouncing. As the phase <^tiit changes fast during the era of 
classical contraction (see (j3.18j) ). it is reasonable to assume that this phase is random 
and uniformly distributed in [0, 2ir) mod 2ir. Under this assumption, we can compute 
the probability of bouncing within the approximation of instantaneous tilt. Namely, 
a bouncing occurs if R vanishes for some time t > itnt- According to (|3.22p . this is 
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equivalent to the condition that there is R € [0, i?tnt) such that 
R 2 



r-, = -UT.max cos 6(R) - p sin(^ tilt ) sin 9{R) (3.23) 

V A z + m z R z 

with #(i?) as given by (|3.21|) . This condition can easily by evaluated numerically. 
Moreover, we obtain the sufficient condition for bouncing that the right side of (13.230 
is negative at R = 0, i.e. 

sm^tiit > „ , (3.24) 

Rtnt P 

where we used that according to (13.2ip . tan#(0) = i?tot/-Rqu- If the right side of (13.241) 
is smaller than one, we can estimate the probability p of bouncing by 

^ 1 ( Rqu |^1, max | 

p > — arccos 1 



7T V -Rtilt P 

Here wi, m ax and p can be rewritten using (13. 3ft or f)3.9|) and (13.15p . 

4. Existence of Bouncing Solutions 

In the previous sections we presented numerical solutions together with approximate 
analytical arguments which indicated that quantum effects can lead to a bouncing of 
the scale function, thus preventing the formation of a space-time singularity. We shall 
now prove that such bouncing solutions exist. The basic strategy is to begin at t = 
at a minimum of the scale function, 

R(0) = , R(0) = , jfj(O) > , (4.1) 

which we choose to be microscopic in the sense that Ri <C Rmt- This solution describes 
an expanding universe. We show that the solution exists all the way to the classical 
era where it goes over to a Friedmann dust solution. 

In order to describe the solution and the initial conditions near the minimum (I4.ip . 
it is useful to rescale the solution and the parameters as follows, 

m ->■ — , t^\, R(t) -> eRit/e 2 ) , A -)• A, w(t) ->■ wit/e 2 ) (4.2) 

EE 1 

(note that this scaling leaves A as well as the length of the Bloch vector unchanged, so 
that the normalization condition (|1 .5|) is preserved). The corresponding rescaled Bloch 
equations are given by 

R\ + e 2 = ~ VA 2 + m 2 R 2 {w e ) u 

(4.3) 



£ 1 



where for clarity, we denoted the solutions by a subscript e. In the limit e \ 0, these 
equations go over to the so-called microscopic limit equations 



I 2 = --^ \A 2 + m 2 R 2 wi , A,, /l 



R = -t^t VA 2 +m 2 R 2 wi , u; = -— ^e 2 Attl. (4.4) 

A z + m z R z 

We remark that that these equations are very similar to the equations (|3.16p used 
in the approximation of instantaneous tilt. The only difference is the additional sum- 
mand +1 in the i?-equation (|3.16p . Also, one should keep in mind that the microscopic 
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limit equations involve the rescaling (|4.2p . whereas the equations (|3.16p merely are an 
approximation of the original system (jl.3p . 



4.1. Solution of the Microscopic Limit Equations. We now solve the approxi- 
mate limit equations for the initial values (|4.ip . First, one sees from the -R-equation 
in (|4.4p that wi(0) = 0. Next, the ^-equation in (|4.4p yields that w 2 is constant. 
Parametrizing w\ and w 3 similar to (|3.19p by an angle 0, we thus obtain 

'wi(t)\ _ ( sin< 

Exactly as described after (|3.19p . the u>-equation can now be solved explicitly, 



w 2 (t) = w 2 (0) 



103(0) with 9\ t=0 = . 



(4.5) 



6{R) = arctan 



\Ri~ 




" R ' 


— arctan 


Rqu 


Rqu 





(4.6) 



In particular, one sees that 0(R) < for all R > Rj, so that wi(t) has the opposite 
sign as 2^3(0). Thus in order to get consistency with the -R-equation in (|4.4p . we must 
choose 

w 3 (0) > . (4.7) 

Then our initial conditions indeed describe a minimum of the scale function (|4.ip . 
More specifically, the i?-equation becomes 

VA 2 + m 2 R 2 



R z 



which can be solved by separation, 



R 2 



sm(@(R)) w 3 (0) 



R dR 

Rr R 



R 



Ri 



VX 2 + m 2 R 2 



sin 



{G(R)) w 3 (0) 



dR. 



This equation is useful for analyzing the behavior as R 
VA 2 + m 2 R 2 



00. Namely, the expansion 



R 2 



shows that 



sin(0GR)) w 3 (0) 



t(R) 



±— 1^3(0) sin arctan 
R 



Ri 
R, 



qu 



+ 



1 

R 2 



Rl 



+ o(Vr) 



m 



(4- 



In particular, our solutions of the microscopic limit equations exist for all times. 

4.2. Construction of Bouncing Solutions. We now turn attention to the rescaled 
equations (|4.3p (for fixed \,m> 0). The continuous dependence of solutions of ODEs 
on the coefficients immediately gives us the following result. 

Lemma 4.1. Let (Ro,wo) be a solution of the microscopic limit equations (|4.4p with 
initial conditions (14. ip and (14. 5p . Then for every T > and S > there is e > and a 
family of solutions (R £ ,w £ ) of the rescaled equations (|4.3p such that for allt £ [—T,T], 

\R £ (t) - R (t)\ < 5 , \w £ (t) - w(t)\ < 5 . 

As we saw in Section STTJ the solution Ro(t) exists for all times, is strictly monotone 
on the time intervals (— 00, 0] and [0, 00) and diverges in the limit t — > 00. The 
function R £ , however, will not be strictly monotone for all positive times. This is 



obvious from the fact that the right side of the -R-equation in (|4.3p tends to zero in 
the limit R — > 00. Thus there must be a minimal time T max (e, Rj) when R £ vanishes. 
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Then clearly R £ does not change signs on the interval [0,T max ]. Moreover, possibly 
by further decreasing e and choosing 5 sufficiently small, we can arrange in view of 
Lemma I4TT1 that T ma _ x (e,Rj) > T. Next, as the angle 9q(R) is negative and decreasing 
according to (|4.6p . we can choose to > 0, k £ (0, 1) and a corresponding tq := Ro(to) 
such that (wq)i(R) < —3k for all R > rp. Then for any T > to, Lemma 14.11 ensures 
that by choosing 6 > sufficiently small, we can arrange that (w E )i(T) < —2k. The 
following lemma shows that by increasing T and further decreasing 5, we can even 
arrange that 



(w £ )i(t) <-k<0 for all t G [T, Jmax] and all sufficiently small e > 



(4.9) 



Lemma 4.2. Let (R £ ,w £ ) be a solution of the rescaled Block equations (j4.3|) defined 
on a time interval [To,Ti]. Assume that R e does not change signs on this interval. 
Then 



arcsm 



N 



arcsm 



K)i(7o) 



N 



< | arctan[i? e (Ti)/i? qu ] — aictan[R £ (To) / R qu ] 
(where N = \w\ again denotes the number of particles). 



(4.10) 



Proof. For ease in notation, we omit the subscript e. The first component of the 
w-equation in (|4.3|) reads 

Ami? 



A 2 + m 2 R 2 



w 3 



Employing the inequality u>3 < \J N 2 — w\ = Ny 1 — (wi/N) 2 , we obtain the estimate 



Wl 



N 



< 



XmR 



A 2 + m 2 R 2 



Vl - {wi/Nf 



and thus 



arcsin(wi(t)/AT) 



< 



arctan[R(t) / R, 



quj 





d 




di 


gn. 





avctan[R(t) / R, 



quj 



where in the last step we used that R has a fixed sign. □ 

The next lemma shows that, possibly by further increasing T and decreasing 5, we can 
arrange that this T max is indeed a local maximum of R £ . 



Lemma 4.3. Assume that (14. 9j) holds and that there is t > T such that 



R £ (t) > 



and 



R £ (t) = 0. 



(4.11) 



Then R £ (t) / 0. 



Proof. We only consider the original ED equations, because then the result is imme- 
diately obtained by rescaling according to (|4.2p . Assume that that statement of the 
lemma is false. Then there is t > T such that P~9|) and (|4.1ip hold but R(t) = 0. 
Combining (|2.5p with (|2.6p and using (|2.7p . we find 



2i?,R = -2(i? 2 + l) + ^(|a| 2 -|/3| i 



(4.12) 
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Again using (|2.5|) . we conclude that at t, 

m 2A _ 

- (|a| 2 - |/3| 2 ) = 2 and Re(a0) = 1 • 

Using (I2.9p . we obtain for the Bloch vector i; that 

2i? R 2 
V3 = — and vi = — — . 

to A 

Computing the corresponding vector w = Uv with [/ as given by (12.101) , we find that 

1 R 2 R R 2 + 2R 2 

Wi = = — and W3 — 



It follows that 



R 2 + R 2 U m JR 2 + R 2 qu 



\wi\ ^ \wi\ R A ^ A i? qu 



mm) J-maxJ 
max 



- |«j 3 j mR(R 2 + 4R 2 U ) ~ mR R ' 
in contradiction to our assumptions (|4.9p and (|4.1ip . □ 

Arguing similarly for negative times, we find that the scale function also has a maxi- 
mum at some time T m i n < —T. We have thus proved the following result. 

Lemma 4.4. For every sufficiently large T and every sufficiently small e > 0, there is a 
solution (R E ,w £ ) of the rescaled Bloch equations (|4,3p defined on the interval [T m i n , T t 
with T m ; n < — T and T max > T with the following properties: At times t = T x 
and t = T m [ Q , the scale function has a local maximum, 

R £ (t) = 0, R £ (t)<0. 

Moreover, the scale function R e is monotone decreasing on [2~min> 0] and monotone 
increasing on [0,T max ]. 

Now we can complete the proof of Theorem 11.11 After performing the scaling (]4.2p 
backwards, we obtain a solution (R, w) of the original equations (ll.3p . By decreasing e 
or increasing T, we can make the scale function and the total time of the contraction 
and expansion cycle as large as we want. We finally note that for fixed A, our construc- 
tion involves the free parameters Ri, w% > as well as any sufficiently small e > 0, 
giving rise to a continuous three-parameter family of solutions. This proves Theo- 
rem [LTJ 

We conclude this section by analyzing the solutions (R £ (t),w £ (t)) of Lemma [44 
outside the time interval [—T,T] in more detail. In particular, we want to show that 
the solutions go over to the Priedmann dust solution as |i| gets large. 

3 

Proposition 4.5. The function Ri is concave if 



e 

Proof. Again, we only consider the unsealed ED equations, because then the result 
follows by rescaling according to (|4.2p . Using (|2.5p . we rewrite (|4.12p as 

2A — 

2RR = -(R 2 + 1) + — Re(a(3) . (4.13) 
R z 

It follows that 

^Rl = I R-k (2RR + R^=^R-h(-l + ^ Re(a^) 
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According to the definition of the Bloch vectors (|2.9p and our normalization conven- 
tion (fT3j) . 

\a\ 2 + |/3| 2 = \v\ = \w\ = N . 

We conclude that 

d 2 3 3 „ i / XN\ 

P2 < - R~2 ( — 1 -I I 

showing that R? is indeed concave if R > XN. □ 

Clearly, this statement implies in particular that R £ itself becomes concave. 

Now we can estimate the time T max (e, Rj) and the corresponding value of R £ from 
above and below. 

Proposition 4.6. There are constants ci, 02,03,04 > such that for all sufficiently 
small e > 0, 

^\<T ia Ue,Ri)<% (4.14) 

J < iimax := i?(T max (e, Ri)) < ^ . (4.15) 

Proof. We choose T so large that the assumptions of Lemma 14.31 are satisfied. Then 
according to (|4.9p and (|1.5p . 

— N < (w £ )i(t) < — k for alU > T (4.16) 

and all sufficiently small e. Using these inequalities in the i?-equation in (|3.16j) and 
setting i?(T max ) = 0, we obtain (|4.15p . 

In order to derive (|4.14p . we note that obviously R e (t) < R(t), where R is the 
solution of the initial value problem 



R = V\ 2 + m 2 R 2 , R(T) = R £ (T) . 
R 2 

By explicit integration we find that R < ct 2 ^ for some c > 0. It follows that 
cT max (e,^/) 2/3 > R(T max (£,Ri)) > ReiT^ie,^)) > c 3 /e 2 , 

giving the lower bound in (|4.14p . The upper bound follows similarly from the inequal- 
ity R £ (t) < R(t) where now R is the solution the initial value problem 

:2 



R = Vx 2 + m 2 R 2 , R(T) = R £ {T). □ 
R 2 

Note that the scalings in Proposition 14.61 coincide with those of the Friedmann dust so- 
lution. This quantifies that for large times, our ED universe indeed behaves classically. 
Our constructions are illustrated in Figure [2j 

4.3. The Energy Conditions. Let us analyze whether our bouncing solutions satisfy 
the energy conditions as described in |13t Section 4.3]. We first bring the energy 
conditions into a form suitable for our analysis. Recall that the weak energy condition 
demands that T a i,W a W b > for any timelike vector W G T p M. Likewise, the dominant 
energy condition states that T ab W a W b > and that T ab W a is a non-spacelike vector. 
Finally, the strong energy condition states that R a t,W a W b > 0. Due to homogeneity, 
it here suffices to evaluate these conditions at a space-time point p with radius r = 0. 
Moreover, due to isotropy, we can assume that the spatial part of W points into the 
radial direction, i.e. 

W = rd t + pd r with < g ab W a W b = t 2 - R 2 p 2 . 
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Figure 2. The bouncing solution R £ . 



Then 



T ab W a W b = T °r 2 - T r r R 2 p 2 

g ac (T b a W h ) (T d c W d ) = r 2 (T ) 2 - R 2 p 2 {T r T f 

whereas the components of the Ricci tensor are computed by 

T = r a = r o + 3r r 



Rt = 8ttk 2T 



Rq 



1 



T 51 



4ttk (r ° 



2 

337) 

R r r = attk (-t ° - t;) . 

Using these formulas, the energy conditions become 
weak energy condition: 
dominant energy condition: 
strong energy condition: Tq > max (317,77). 

Using (|2.9p . we can express (|2.4p in terms of the Bloch vector v, 



Tq > max 
T ° > 1 77 1 



0,77) 



rpO 



A 

m v s - — vi 



1 r 



Computing the vector w = Uv as in the proof of Lemma I4.3| we obtain 



rpO 



m 



-^4 \J R2 + R \n Wl 



3mRA Jr 2 + r 2 



3R 3 



qu 



R 2 + R 2 U 



w 3 . 



Rescaling by (|4.2p only gives an overall factor of e 

m 



4 T 



e 4 TI 



Ri 



A- 1 



3mR e JR 2 + R 2 



(W £ )l + 



qu 



3R I < IR 2 + R 2 



qu 



(4.17) 



(4.18) 



(4.19) 
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We first show that the energy conditions are satisfied for large R £ , as was to be 
expected in view of the Friedmann dust limit. 

Proposition 4.7. Assume that f|4.9[) holds and that 

Re > • (4.20) 

K 

Then the weak, dominant and strong energy conditions are satisfied. 

Proof. We again only consider the unsealed equations. Using (|1.5|) and (|4.9p . we can 
estimate (j4.18j) by 

Hence 

where in the last step we applied (|4.20p . □ 

A more interesting question is whether the energy conditions are satisfied in the 
bouncing regime. We first prove our results and discuss them afterwards. 

Proposition 4.8. Suppose that a solution bounces at time t, i.e. 

R £ {t) = and R £ {t) > . 

Then the strong energy condition is violated at t. 

Proof. We again consider the unsealed equations. Comparing (|4.18p with (|1.3p . one 
sees that 

T„° = ^ + l) = i. 
Moreover, combining (14.13P with (12.4p . we find 

2A — 

< 2RR = -1 + — Re(a/3) = -1 + 3R 2 Tl 
R z 

and thus 

in contradiction to the strong energy condition in (|4.17p , □ 



Proposition 4.9. Let (R £ ,w £ ) be a the solutions of Lemma \4-4\ Then for sufficiently 
small e > 0, the weak, the dominant and the strong energy conditions are violated at 
the bouncing time. 

Proof. Suppose that R £ (t) = 0. Solving the i? £ -equation in (|4.3p for (w e )i and substi- 
tuting the result into (|4.19p . we find that 

1 



^0 _ -2 



Rl 



rpr _ _ -2 ^qu I -4 _^ \ 

m^Ri + R\ v ) + sri^^rT 3 - 

For small e, the last summand of TJT dominates, which in view of (|4.7p is positive. □ 
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The statement of Proposition 14.81 rules out any bouncing solutions which satisfy the 
strong energy condition. This also seems necessary for consistency, because otherwise 
such bouncing solutions would violate the Hawking-Penrose singularity theorems (see 
in particular |X3|, Theorem 2 in Section 8.2]). On the other hand, the singularity 
theorems do not rule that bouncing solutions might satisfy the dominant or weak 
energy conditions. Proposition 14.91 shows that the methods of Section [4.21 do not give 
rise to a procedure to construct such bouncing solutions. However, it is still possible 
that such solutions do exist in a certain parameter range where e is bounded away 
from zero. Indeed, by choosing at the bouncing time positive but sufficiently small, 
one can arrange in view of (|4.18p that the dominant and weak energy conditions are 
satisfied near the bouncing point. A numerical study showed that for such initial 
conditions, a typical solution does not enter the classical dust region, but instead the 
scale function has a maximum after a short time and then decreases again, leading 
to a singularity. The question whether there are special initial conditions where the 
solution does enter the classical dust region seems a hard global problem which we 
cannot enter here. 

5. Existence of Time- Periodic Solutions 

For the construction of time-periodic solutions we shall make use of the following 
symmetry argument. 

Lemma 5.1. Let (R e ,w e ) be a solution of the rescaled Block equations (|4,3p defined 
on a time interval [Tq,T\]. Assume that for some time T in this interval, 

R £ (T)=0 and (w £ ) 2 (T) = . (5.1) 

Then for all t with T ± t € [To, T±], 

R £ (T + t) = R £ (T - t) and w £ {T + t) = P w e {T - t) , 

where the operator P flips the sign of the second component, 





(wi\ 




( m \ 




U>2 


H 


-w 2 








{ w 3 J 



Proof. One easily verifies that both (R £ {T + t),w(T + t)) and {R £ (T - t),Pw{T - t)) 
are solutions of the rescaled Bloch equations (I4.3P (note that the corresponding R- 
equations are the same, whereas in the w-equations the second term in the brackets 
on the right has opposite signs; the signs of the components of w are chosen such the 
i/J-equations hold in both cases). According to (15. lh . these solutions have the same 
initial values at t = 0. Therefore, we would like to apply the uniqueness theorem for 
ODEs to conclude that the solutions coincide. However, this is not possible directly 
because the i?-equation in (|4.3p is singular if R = 0. In order to bypass this problem, 
we consider the corresponding vectors (R, R, w) as solutions of the system of ODEs 
consisting of the w-equation in (|4.3p . the equation j^R = R and the second order 
equation (|2.6p (also rescaled by (|4.2p ). Then the singularity at R = has disappeared, 
making it possible to apply standard uniqueness results. □ 

Proof of Theorem M.SX For given A and Ri, we let (R e ,w £ ) be a family of solutions of 
the rescaled Bloch equations (|4.3p with initial conditions 

R £ (0) = R I , R £ (0)=0 and (^ e ) 2 (0) = , (w £ ) 3 (0) > . 
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Following the constructions of Section H~2l we know that for every sufficiently small e > 
0, we have a bouncing solution in the sense of Lemma 14.41 Moreover, according to 
Lemma 15.11 we know that 

R £ {-t) = R £ {t) and w £ (-t) = Pw £ {t) . (5.2) 

Hence there is T max > where the scale function has a maximum. If (w e )2(T) vanishes, 
we see from (|5.2|) that our solution (R £ ,w £ ) has the same values at t = ±T max . Hence 
the uniqueness theorem for ODEs yields that the solution can be extended to all times 
and is 2T max -periodic. Thus it remains to show that for any Ri, we can choose e > 
such that (w £ )2(T) = 0. 

Parametrizing w £ similar to (|3.14p by w £ = {{w £ )\,p £ cos 4> £ , p £ sin cj) £ ) , this condition 
is equivalent to demanding that 4> £ (T max (e)) £ tt/2 + irZ. In order to arrange this 
condition, we recall that R £ (t) is continuous in e for each fixed t. Since R £ is concave 
around its maximum value for e small enough, we see that R mgJ i(£) is continuous 
in e, hence so is T max (e). As a consequence, w £ (T max ) and thus also £ (T max (e)) is 
continuous in e. Next, combining the bounds (14,14p and (I4.16p . one sees from the 
w-equation that </> £ (T max (e)) tends to infinity as e \ 0. Thus from the intermediate 
value theorem, there is e such that </> e (T max (e)) G tt/2 + irZ. □ 

Note that the last proof even shows that for every fixed Ri, the admissible e form 
an infinite series (e n )neN with lirm^oo e n = 0. We also remark that additional time- 
periodic solutions could be constructed by choosing (iu e )2(0) / and adjusting Ri 
and e so as to arrange that {w £ )2{T m \ n ) = = (w £ )2(T mgx ). Then in general R £ (T m i n ) ^ 
Re(T m ax), but by applying Lemma ISTTI one sees that the solution is still time-periodic 
with the doubled period 2(T max — T m ; n ). 

6. Solutions in an Open and Flat Universe 

So far, we only analyzed the closed case. We now outline which of our methods 
and results carry over to the the open and flat cases. We first point out that the 
parameter k which distinguishes the different cases only enters the 72-equation in (jl.3p 
via the summand +k. In the closed case, the summand +1 is essential for the universe 
to form an equator in the classical dust limit. Thus in the open and flat cases, there is 
no mechanism which contracts the universe after entering the classical era. However, 
in our analysis of the quantum era, the summand +/c in the /^-equation was never 
used. Indeed, it even dropped out of the equations in the microscopic limit (|4.4p . In 
particular, our construction of bouncing solutions applies just as well in the open and 
flat cases. The only difference is that after entering the classical era, the universe will 
keep expanding forever. This is illustrated in Figure [3] in a typical example. Clearly, 
no time-periodic bouncing solutions exist. 

A specific difference in the open case is that the component w\ can be positive 
(see the -R-equation in (jl.3p ). In view of (I4.18p . this gives rise to bouncing solutions 
for which the energy conditions are violated even for large times (see Figure H] for a 
typical example). These solutions correspond to the unphysical situation of a negative 
classical energy density which drives an accelerated expansion of the universe. 

Appendix A. The Dirac Equation and its Separation 

For the derivation of the corresponding Dirac operator we follow the methods in [7], 
for an alternative derivation see |18j . For the Dirac matrices in Minkowski space we 
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Figure 3. Example of bouncing solutions in a closed, flat and open 
universe for m = 21.5, A = 3/2, R max = 0.1, i? ti it » 0.0698. 




FIGURE 4. Example of a bouncing solution in the open case with neg- 
ative energy density for m = 21.5, A = 3/2, i? max = lj Rtttt ~ 0.0698. 
For this solution, w\ ~ 0.046 is almost constant. 



always choose the Dirac representation 




where a a , a = 1,2,3, are the Pauli matrices, 




(A.l) 



In curved space-time, the Dirac wave functions are four-component complex functions 
on the manifold, which at every space-time point are endowed with an inner prod- 
uct of signature (2,2). For convenience, we represent this inner product just as 
in Minkowski space as ^ = ^^/°^, thus ^ is the usual adjoint spinor. The first step 
for the construction of the Dirac operator is to choose Dirac matrices G 3 which are 
symmetric with respect to the inner product and satisfy the anti-commutation 
relations 



{G j ,G k } = G j G k + G k G j = 2g jk U 



(A.2) 
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This can be accomplished by taking the following linear combinations of the Dirac 
matrices of Minkowski space, 

m 



G 
G 

G* 



R(t) v 
1 

rR(t) 
1 



cos i? 7 3 + sin i? cos p> 7 1 + sin 1? sin ip 7 2 ) 
(— sin -d 7 3 + cos $ cos p 7 1 + cos 1? sin p 7 2 ) 
(— sin 93 7 1 + cos ip 7 2 ) . 



(A.3) 



rR(t) sin?? 

The spin connection D can then be written as 

Dj = Oj - iEj , 

where Ej, the so-called spin coefficients, are given in terms of the Dirac matrices and 
their first derivatives. Thus the Dirac operator becomes 

V = iGWj = iG j dj + (GPEj) . 

For an orthogonal metric, the combination G 3 Ej takes the simple form (for details 
see [9l Lemma 2.1]) 



GPE~ 



dj ( \f\9\G j ) with g = det gij 



making it unnecessary to compute the spin connection coefficients or even the Christof- 
fel symbols. A short calculation gives 

11 [ at + 2R(t) J + R(t) \-V n J ' 
where the purely spatial operator 2?-^ is given by 



(A.4) 



Vh ia< [dr + ^j^) +icr"i) ll + iar<) r .. 



(A.5) 



and the matrices a a , a £ {x, $,ip}, are the following linear combinations of the Pauli 
matrices, 

® T := f( r ) (cosfi a 3 + sin'!? cosy? a 1 + siwd sin ip <r 2 ) 
* 1 



(7 



sin i? cr 3 + cos d cos 1^ a 1 + cos 1? sin ip a 2 ) 



1 



sin ip a 1 + cos </> <r 2 ) . 



r sin$ 

In order to separate the Dirac equation 

(V - m)ty = , 
we consider an eigenspinor -0a of the spatial Dirac operator, 

and employ the ansatz 



(A.6) 



(A.7) 



(A.« 
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we obtain the coupled system of ODEs (|2.2p for the complex- valued functions a and (3. 
Using that the matrix in (|2.2p is Hermitian, one easily verifies that 

|(l«| 2 + l/3| 2 )=0. 
Thus we can normalize the spinors according to (|2.3p . 

Appendix B. The Energy-Momentum Tensor 

We begin in with the simpler closed case (the flat and open cases will be treated in 
the paragraph following (|B.6p below). Then the operator T>^ = T> S 3 with domain of 
definition C°°(S^) 2 ) is an essentially self-adjoint operator on L 2 (S 3 ) 2 . It has a purely 
discrete spectrum with eigenvalues A = ±|,±|, .. . and corresponding eigenspaces of 
dimension N = A 2 — \ (see for example (9j Appendix A]). For a given eigenvalue A 

we let ifis*~' , . . . , be an orthonormal basis of eigenvectors. Using the ansatz (|A.8[> 
(for fixed functions a and /3), we introduce corresponding one-particle wave func- 
tions ^i, . . . , 'I' at. Taking their wedge product, 

q> ■— q> 1 a • • • A *tv , 

we obtain a state of the fermionic Fock space. This state is a homogeneous and 
isotropic Hartree-Fock state. 

The energy-momentum tensor is introduced as the following expectation value, 

T jk = X - Re&KiGjDk + iG k Dj)*)s , 

where is the scalar product on the Fock space associated to the probability 

integral, and the operator (iGjD k + iG k Dj) acts on the Fock space as a one-particle 
operator O, i.e. 

O^/ = A • • • A ^> N + #i A (C* 2 ) A • • • A ^> N H h f i A ■ • ■ A (OV N ) . 

Using the orthonormality of one-particle wave functions ^j, the above expectation 
value reduces to the sum of the one-particle expectations, 

1 - — 

T ik = - Re (iGjD k + iG k Dj) * a . 

1 a=i 

In order to bring the energy-momentum tensor into a more explicit form, it is 
convenient to introduce to introduce the object 

P(t, x; t', x') = *a(t, Z)W, ■ (B.l) 

a 

Denoting the spectral projectors of P53 by E\ and keeping in mind that the adjoint 
spinor in (jB.ip involves multiplying by the matrix 7 , we can write (IB.ip as 

The energy-momentum tensor can be expressed in terms of P by 

T jk (t, x)= l - Tr c4 {(iGjD k + iG k D 3 )P(t, x; t', ^))\ t , =t . l=3 , (B.3) 

where the trace is taken over 4 x 4-matrices. Now we could proceed by computing 
the spin connection and substituting the resulting formulas into (|B.3p . The following 
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shorter method avoids the computation of the spin coefficients: Collecting the time- 
dependent part in (|A.4|) . we see right away that D t = dt — 3R/(2R). Hence 



t'=t,x' 



and employing the ansatz (|B.2p . we obtain 



T ° = R~ 3 



™(M 2 " l/5| 2 ) - -^-Re(a/3)J Tt C 2E x (x,x) . (B.4) 

Next it is easy to compute the trace of the energy- momentum tensor. Namely, using 
that V = iG j Dj, we get from (lR3|) 

T£ = Tr C 4 {VP(t,x;t',x')} \ v=t ~ t=s = mT ¥ {P(t,x;f ,x')} \ v=t ~ l=s 

= R- 3 m(\a\ 2 - |/3| 2 ) Tr c2 £ A (x,x) . (B.5) 

The quantity Tr C 2 E x (x,x) appearing in (jB.4[) and (|B.5p is a constant. It is computed 
by 

Tr c2 E x (x, x) = Tr c2 E x (x, x) d^x) = tv(E x ) = (a 2 - I) . 



In view of our normalization convention (|2.3j) for a and /3, it seems most convenient 
to use the convention n(S 3 ) = 1 and to set 

Tr c2 £ A (£,£) = l. (B.6) 

Knowing (]B.4p and ()B.5P and using the convention (|B.6p . we can determine all compo- 
nents of the energy-momentum tensor by a symmetry consideration. First, spherical 
symmetry implies that T r ° = T® = T° = 0. Next, using that our system is homogeneous 
in space, we conclude that T£ = = = 0. Thus the energy-momentum tensor 
must be a diagonal matrix. Using once again that our system is homogeneous, we even 
obtain that T£ = Tf = . Hence the energy- momentum tensor has only two indepen- 
dent components, which are determined from the two equations (|B.4p and ()B.5|) . We 
thus obtain the equations (|2.4p . and all other components of the energy- momentum 
tensor vanish. 

In the open and flat cases, the situation is a bit more difficult because the opera- 
tor P% has a continuous spectrum, making it necessary to proceed as follows. After 
choosing a self-adjoint extension of 2?%, we can apply the spectral theorem 



V H = A dE x , (B.7) 

Ja(V H ) 

where dE x is the spectral measure. Using that the spectrum is absolutely continuous, 
we can represent the spectral measure by an integral kernel, i.e. 

dE x = E x dX with (E x ^)(x) = J E x (x, y) if>(jfl dfi(y) . 

This integral kernel is no longer a projector. It is an operator of infinite rank, cor- 
responding to the fact that our system now involves an infinite number of particles. 
Also, the states in the image of E x are no longer normalized, but we must rely on the 
distributional normalization 

ExE^ = 5(X -fi)E x . 
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Describing the fermions our our system by E\ (x, y) , all local quantitites (like the 
energy-momentum tensor, the current, etc.) are well-defined. But the total energy, 
charge, etc. will clearly diverge because the spatial volume is infinite. With this in 
mind, we can work with E\(x, y) just as if it were a projector onto the fermionic states 
of our system. Then the above considerations again apply, and for convenience we 
take (|B.6|) as a normalization convention. This again gives the equations (\2A\i . 
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